First-order phase transition in easy-plane quantum antiferromagnets 
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Quantum phase transitions in Mott insulators do not fit easily into the Landau- Ginzburg- Wilson 
paradigm. A recently proposed alternative to it is the so called deconfined quantum criticality 
scenario, providing a new paradigm for quantum phase transitions. In this context it has recently 
been proposed that a second-order phase transition would occur in a two-dimensional spin 1/2 
quantum antiferromagnet in the deep easy-plane limit. A check of this conjecture is important for 
understanding the phase structure of Mott insulators. To this end we have performed large-scale 
Monte Carlo simulations on an effective gauge theory for this system, including a Berry phase term 
that projects out the S — 1/2 sector. The result is a first-order phase transition, thus contradicting 
the conjecture. 
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The Landau-Ginzburg- Wilson (LGW) theory for 
phase transitions has been an immensely successful 
paradigm of physics for the last 50 years. It is one 
of the cornerstones of statistical and condensed matter 
physics, providing deep insight into phase transitions [l[ . 
The standard example is the well-known paramagnetic- 
ferromagnetic phase transition. Recently, examples of 
phase transitions that do not fit into the LGW paradigm 
have been discussed @, H, 0| • A prominent example is the 
continuous quantum phase transitions from a Neel state 
with conventional antiferromagnetic order into a param- 
agnetic valence-bond solid (VBS) state 0]. In the Neel 
state, an SU(2) symmetry is broken, while in the VBS 
phase translation invariance of the lattice is broken. The 
LGW paradigm does not describe this phase transition 
correctly, since it predicts a first-order phase transition in 
this case. In view of the failure of the LGW paradigm in 
this and other cases, a new scenario has recently been 
proposed 0, 0], introducing the concept of deconfined 
quantum criticality (DQC). This concept applies to sys- 
tems where the order parameter can be viewed as being 
composed by elementary building blocks. For instance, 
in the case of the Neel- VBS transition the spinons are 
the building blocks of the spin field. Similarly to quarks 
in hadrons, the spinons are confined in both the Neel 
and VBS phases. The DQC scenario asserts that the 
spinons are deconfined only at the critical point. This 
claim is based on a subtle destructive quantum interfer- 
ence mechanism between instantons and the Berry phase 

1- 

An attempt to provide proof of evidence for these ideas 
has recently been put forth |2| ■ It involves a deformation 
of the two-dimensional Heisenberg model into an easy- 
plane quantum antiferromagnet. The effective theory of 
a spin 1/2 quantum antiferromagnet is a 0(3) nonlinear 
a model with a staggered Berry phase factor [1] . Such a 
nonlinear a model describes the fluctuations of the ori- 



entation \\j of the order parameter. The easy-plane de- 
formation adds a term proportional to t? z a to the action, 
which explicitly breaks the 0(3) symmetry down to U (1). 
This lower symmetry simplifies considerably the analysis, 
especially when the CP 1 representation rij = z* a a a bZjb 
is used, with \zji\ 2 + \zj2\ 2 = 1 due to the local constraint 
nj = 1. The CP 1 representation naturally introduces a 
local abelian gauge symmetry, since is invariant un- 
der the local gauge transformation zj a — ► e ldai zj a . A deep 
easy-plane deformation forces n 2 j w 0, thus inducing the 
additional local constraint \zji\ 2 w |zj2| • This allows us 
to write Zj a — e i a /^/2. The requirement of local U(l) 
gauge invariance and the deep easy-plane limit naturally 
leads to an effective lattice gauge theory for a quantum 
antiferromagnet proposed in Ref. [fj, which for S = 1/2 
has the lattice Lagrangian 
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Cj = -f3^2cos(A tl 9j a -Aj^)-Kcos(e flv \A u A j x)+irijAj T , 

a=l 

(1) 

where Aj is a compact gauge field which here is doing 
more than just being an auxiliary field, like in the case 
of the CP 1 model. It determines also the Berry phase for 
the above model. The index r corresponds to imaginary 
time and the staggering factor is given by r]j = (— 1) J . 
Note that the present gauge field is a function of the 
spacetime coordinates, in contrast to the usual Berry 
gauge potential appearing in spin models which is 
a functional of the spin field. 

Compactness of the gauge field gives rise to instan- 
ton configurations which are known to spoil the phase 
transition in a corresponding model with only one phase 
field and no Berry phase [8j. In the absence of Berry 
phase and with two phase fields present, on the other 
hand, aphase transition in the 3DXY universality class 
occurs [§, [13] regardless of whether the gauge field is 
compact or not. Since only one gauge field is present 
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and there are two phase fields available, the Higgs mech- 
anism is able to supress only one out of two massless 
modes. The remaining massless mode is charge neutral 
and drives the ZDXY transition [loj ]. 

Recently, it has been argued that the Berry phase, 
which is crucial to describe the phase inside the paramag- 
netic phase @,[3], suppresses the instantons at the critical 
point. Here we will investigate this point by monitoring 
the phase transitions in the model (Q} in the presence 
and absence of a Berry phase. In the former case a phase 
transition is expected in the charged sector, contrasting 
with the transition driven only by the neutral sector in 
absence of the Berry phase. The result will be shown to 
be a first-order phase transition. 

The DQC scenario implies that the critical point is 
governed by an easy-plane system Lagrangian featuring 
a non-compact gauge field, i.e., 
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Ci = - / 3y^cos(A M 6'j, 



^)+2<AxA,f 
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This model with unequal bare phase stiffnesses has been 
studied in great detail (lQ |. It features two distinct 
second-order phase transitions, one belonging to the 
3DXY universality class and another one correspond- 
ing to the so-called inverted 3DXY transition In 
the limit where the bare phase stiffnesses are equal, clear 
signals of non-SDXY behavior are seen 

[13, EE El- In 

Ref. HI strong indications of a first-order phase tran- 
sition in a loop-gas representation of the non-compact 
model Eq. ([2]), were found. We will consider both Eqs. 
(P) and © in detailed Monte Carlo (MC) simulations. 

For performing MC simulations on the model with a 
Berry phase term, it is convenient to introduce a dual 
representation of the model Eq. (TTJ). In such a repre- 
sentation the action is real, with a Lagrangian given by 



A large number of sweeps is required in order to get ad- 
equate statistics in the histograms (see below) for the 
largest system sizes. Firstly, we have computed the 
second moment of the action M-i = ((S — (S)) 2 ) for 
the model with and without a Berry phase term, where 
S = J2i^i- Secondly, we have focused on a number 
of quantities that provide information on the character 
of the phase transition associated with the specific heat 
anomaly. The first of these quantities is the third mo- 
ment of the action, M3 = ((S— (S)) 3 }. At a second-order 
phase transition this quantity should scale as follows. 
The peak-to-peak height scales as i( 1+Q )/ !y whereas the 
width between the peaks scales as L~ x l v 14]. At a first- 
order phase transition, these quantities scale as L 6 and 
L -3 , respectively [13]. We also study the probability dis- 
tribution P(S, L) of the action S for various system sizes. 
At a first-order phase transition, P(S, L) will exhibit a 
double-peak structure associated with the two coexisting 
phases. 




FIG. 1: (Color online) Specific heat M2 of Eq. Q for various 
system sizes. Panel (a): Without Berry phase term. The 
peak develops into a singularity of the 3DXY type. Panel (b): 
With Berry phase term. The peak develops into a 5-function 
singularity with a peak scaling as L A , consistent with a first- 
order transition. Note the symmetry and asymmetry of the 
peaks in the right and left panels, respectively. This is to be 
expected, since the peaks in the right panel originate with the 
superposition of a 3DXY peak and an inverted 3DXY peak. 



.(2) 



-f l + As l ) 2 - (3) 



Here, h( a ) are integer-valued dual gauge fields, and 
E^lSy = S^rji- Note that we would obtain Eq. ^ 
both for Eqs. JT) and ©, with /* = for Eq. ©. For 
Eq. {1} with compact A^, Sj is integer- valued. For Eq. 
([2|) with a non-compact A^, Sj is real- valued. There- 
fore in the former case Sj can be gauged away since the 
h( a '-fields are integer-valued. We have chosen a gauge 
where = 0. The MC computations were performed 
using Eqs. © and ©. For both Eqs. © and ©, wc 
have used k = (3. We have used the standard Metropolis 
algorithm with periodic boundary conditions on a cubic 
lattice of size L x L x L. For Eq. ((3]) we have used 
L = 4,8,12,16,20,24,32,36,48,60,64,72,80,96,120, 
while for Eq. © we have used L = 48, 64, 80, 96, 112, 120. 



The specific heat M2 is shown in Fig. [TJ Panel (a) 
shows the anomaly for the model Eq. © with no Berry- 
phase term, i.e., A x f = (0, 0, 77) = 0. The anomaly has 
the characteristic asymmetric shape of the SDXY model. 
In this case, there are no Berry phases to suppress the 
instantons of the compact gauge-field A at the critical 
point. Hence, the charged sector does not feature critical 
fluctuations that can interfere with those of the neutral 
sector. When the Berry phase field f is included, the 
specific heat is notably more symmetric and the anomaly 
develops into a (5-function peak, consistent with a first- 
order phase transition. This is shown in panel (b). 

To investigate more precisely the character of the phase 
transition when a Berry phase term is present, we have 
performed finite-size scaling (FSS) of the third moment 
of the action, M3 14|. The results are shown in Fig. 



panel (a). It is seen that for small and intermediate 
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FIG. 2: (Color online) Scaling of the height [panel (a)] and 
width [panel (b)] of M3 of the action in Eqs. ^ and (3J. 
The lines in panel (a) represent L 1 ' 43 and L 6 . The former 
is the 3DXY result. The lines in panel (b) represent L -1 ' 49 
and L~ 3 . The former is the 3DXY result. For large system 
sizes, the height and width scale in manner consistent with 
a first-order phase transition. Also shown are results for Eq. 
© with no Berry phase term f = (green symbols). These 
results follow the 3DXY scaling lines. The red symbols are 
the results for Eq. (J3J) while the blue symbols are results for 
Eq. ©. 



system sizes, the height increases with L in a manner 
which might appear consistent with that of a second- 
order phase transition. However, the quality of the scal- 
ing is not satisfactory, since a clear curvature in the scal- 
ing plots is seen (red data points). As system sizes in- 
crease we see a gradual increase in the apparent value of 
(1 + cx)/v, until for large system sizes, we clearly have 
A/a r~j L e , consistent with a first-order phase transition 



Panel (b) of Fig. [5] shows the scaling of the width of 
M3. Again, the line with the smallest negative slope is 
the line one would obtain for the 3DXY model, while 
the line with the most negative slope is ~ L~ 3 , charac- 
teristic of a first-order phase transition. Again we obtain 
apparent scaling, with a crossover regime at intermediate 
length scales into a regime where the width scales as it 
would in a first-order phase transition [15J. The results 
of Fig. [2] provide further support to the notion that the 
phase transition in the model with a compact gauge field 
and a Berry phase term is a first-order phase transition. 

To investigate this further, we have computed the 
probability distribution P(S, L) for various system sizes. 
The results are shown in Fig. [3] Panel (a) shows re- 
sults for Eq. (JTJ) in the representation Eq. ([3]). Panel 
(b) shows results for Eq. The Ferrenberg-Swendsen 
algorithm has been used to reweight the histograms [lj| • 



For L < 48, we essentially have not been able to resolve 
a double peak structure at all, showing that the phase 
transitions in the models Eqs. dXJ) and @ are weakly 
first-order. We have located the transition temperature 
from the peak structures in the specific heat and M3, and 
performed long simulations at this temperature for each 
L. For the largest systems, L = 96, 120, up to 120 • 10 6 
sweeps over the lattice were done. A clear double-peak 
structure in P(S, L) is seen to develop for system sizes 
L > 60. The fact that such large system sizes are re- 
quired to bring out the double-peak structure, implies 
that this phase transition is weakly first-order. 
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FIG. 3: (Color online) Histograms for the probability distri- 
bution P(S, L) as a function of S/L 3 for various system sizes 
L. (a): results for Eq. (fl| in the representation Eq. J3J. (b): 
results Eq. @. A double peak structure develops with the 
latent heat per unit volume approaching a finite constant as 
L is increased. This is a hallmark of a first-order transition. 
For the largest systems, up to 120 • 10 6 sweeps over the lattice 
were performed. A total of approximately 500000 CPU hours 
were used to obtain these results. 



We also perform FSS of the height of the peak be- 
tween the two degenerate minima in the free energy 
— ln[P(S, L)]. This height should scale as L 2 in a first- 
order phase transition, since it represents the energy of 
an area which separates two coexisting phases |17| . The 
results are shown in panel (a) of Fig. |4j For large enough 
systems, the height clearly approaches the dotted line 
~ L 2 , as in a first-order transition. This is corroborated 
by extracting the latent heat per unit volume in the tran- 
sition, shown in the lower panel of Fig. [¥] It approaches 
a nonzero constant as L is increased, as it should in a 
first-order phase transition. 

Further insight into the nature of the first-order phase 
transition can be obtained by means of the renormaliza- 
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transition in a closely related model has also been found. 

In summary, our large scale MC simulations of the 
deep easy-plane quantum antiferromagnet confirm the 
instanton-Berry phase suppression mechanism proposed 
in Ref. Therefore, the spinons are indeed deconfined 
at the phase transition. However, the phase transition in 
this case is first-order, which contradicts the DQC picture 
for this model, where a second-order phase transition has 
been predicted. 
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(NANOMAT), and 167498/ V30 (STORFORSK), and 
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FIG. 4: (Color online) Panel (a) shows the scaling of 
the height AF of the peak between the two minima in 
- In P(S,L) both for Eqs. © (red curve) and © (blue 
curve). Dotted line is the line ~ L 2 . The height scales as 
AF ~ L A ' X . This is a hallmark of a first-order transition. 
Panel (b) shows latent heat per unit volume AS/V as a func- 
tion of L. The upper (blue) curve is from Eq. ([2}, the lower 
(red) curve is from Eq. ©. AS/V approaches a nonzero 
value as L is increased. 



tion group (RG). In the field theory Lagrangian the in- 
teraction of an easy-plane system reads C lBil — u {\zi\ 2 + 
\z 2 \ 2 ) 2 /2 + v \ Zl \ 2 \z 2 \ 2 = U0 (N 4 + M 4 )/2 + W0 N 2 M 2 , 
where wq = uq + vq. Consider a generalized situation 
where the complex fields have each N/2 components. 
The renormalized dimensionless couplings in d = 4 — e 
dimensions are g — ufi~ £ , h — uifi~ E , and /, where / is 
the dimensionless gauge coupling and [i is an arbitrary 
mass scale. The j3 functions at one-loop order are [19| 
[3 g = -eg - Ggf + (N + 8)g 2 /2 + 2Nh 2 + 6/ 2 , [3 h = 
-eh - 6hf + 3 (AT + 2)gh + 6/ 2 , and [3 f = -ef + Nf 2 /3. 
Nontrivial fixed points with / = 3e/N and h < are 
found for N > 300, while in the deep easy-plane limit 
h = no fixed points with f = 3e/N are found for all 
values of N. In a Ginzburg-Landau (GL) theory of super- 
conductors, the existence of a critical value of N above 
which nontrivial fixed points are found actually reflects 
the strong-coupling behavior at much lower values of N. 
It turns out that the phase transition for one complex 
order parameter is second-order in the type II regime 
111 ], while a first-order transition occurs in the type I 



regime [2fJ, |2l|. Inspired by the GL case, we interpret 



the complete absence of a critical value of N for h = as 
a clear signature of a first-order phase transition in the 
deep easy-plane regime. This is a further confirmation 
of our MC results; see also Ref. [12|], where a first-order 
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